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O ' Abstract 

o : 

' Equal-arm interferometric detectors of gravitational radiation allow phase measurements many 

orders of magnitude below the intrinsic phase stability of the laser injecting light into their arms. 
This is because the noise in the laser light is common to both arms, experiencing exactly the 
same delay, and thus cancels when it is differenced at the photo detector. In this situation, much 

■ lower level secondary noises then set overall performance. If, however, the two arms have different 

m ! 

, lengths (as will necessarily be the case with space-borne interferometers), the laser noise experiences 

' different delays in the two arms and will hence not directly cancel at the detector. In order to 

■ solve this problem, a technique involving heterodyne interferometry with unequal arm lengths and 

Q^i independent phase-difference readouts has been proposed. It relies on properly time-shifting and 

5-H i 

linearly combining independent Doppler measurements, and for this reason it has been called Time- 
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Delay Interferometry (or TDI). This article provides an overview of the theory and mathematical 
foundations of TDI as it will be implemented by the forthcoming space-based interferometers such 
as the Laser Interferometer Space Antenna (LISA) mission. We have purposely left out from 
this first version of our "Living Review" article on TDI all the results of more practical and 
experimental nature, as well as all the aspects of TDI that the data analysts will need to account 
for when analyzing the LISA TDI data combinations. Our forthcoming "second edition" of this 
review paper will include these topics. 
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I. INTRODUCTION 



Breakthroughs in modern technology have made possible the construction of extremely 
large interferometers both on ground and in space for the detection and observation of grav- 
itational waves (GW). Several ground based detectors are being constructed or are already 
operational around the globe. These are the LIGO and VIRGO interferometers, which have 
arm-lengths of 4 km and 3 km respectively, and the GEO and TAMA interferometers with 
arm-lengths of 600 m and 300 m respectively. These detectors will operate in the high fre- 
quency range of GW of ~ 1 Hz to a few kHz. A natural hmit occurs on decreasing the lower 
frequency cut-off of 10 Hz because it is not practical to increase the arm-lengths on ground 
and also because of the gravity gradient noise which is difficult to eliminate below 10 Hz. 
However, VIRGO and future detectors such as the advanced LIGO, the proposed LCGT in 
Japan and the large European detector plan to go to substantially below 10 Hz. Thus, in 
any case, the ground based interferometers will not be sensitive below the limiting frequency 
of 1 Hz. But on the other hand, there exist in the cosmos, interesting astrophysical GW 
sources which emit GW below this frequency such as the galactic binaries, massive and 
super-massive black-hole binaries etc. If we wish to observe these sources, we need to go 
to lower frequencies. The solution is to build an interferometer in space, where such noises 
will be absent and allow the detection of GW in the low frequency regime. LISA - Laser 
Interferometric Space Antenna - is a proposed mission which will use coherent laser beams 
exchanged between three identical spacecraft forming a giant (almost) equilateral triangle of 
side 5 X 10^ kilometers to observe and detect low frequency cosmic GW. The ground based 
detectors and LISA complement each other in the observation of GW in an essential way, 
analogous to the way optical, radio, X-ray, 7-ray etc., observations do for the electromag- 
netic spectrum. As these detectors begin to operate, a new era of gravitational astronomy 
is on the horizon and a radically different view of the universe is expected to emerge. 

The astrophysical sources that LISA could observe include galactic binaries, extra-galactic 
super-massive black-hole binaries and coalescences, and stochastic GW background from the 
early universe. Coalescing binaries are one of the important sources in the LISA frequency 
band. These include galactic and extra galactic stellar mass binaries, and massive and super 
massive black-hole binaries. The frequency of the GW emitted by such a system is twice 
its orbital frequency. Population synthesis studies indicate a large number of stellar mass 
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binaries in the frequency range below 2-3 mHz [l|, In the lower frequency range (< 1 
mHz) there are a large number of such sources in each of the frequency bins. Since GW 
detectors are omni-directional, it is impossible to resolve an individual source. These sources 
effectively form a stochastic GW background referred to as binary confusion noise. 

Massive black-hole binaries are interesting both from the astrophysical and theoretical 
points of view. Coalescences of massive blackholes from different galaxies after their merger 
during growth of the present galaxies would provide unique new information on galaxy 
formation. Coalescence of binaries involving intermediate mass blackholes could help un- 
derstand the formation and growth of massive blackholes. The super massive black hole 
binaries are strong emitters of GW and these spectacular events can be detectable beyond 
red-shift of 2; = 1. These systems would help to determine the cosmological parameters in- 
dependently. And, just as the cosmic microwave background is left over from the Big Bang, 
so too should there be a background of gravitational waves. Unlike electromagnetic waves, 
gravitational waves do not interact with matter after a few Planck times after the Big Bang, 
so they do not thermalize. Their spectrum today, therefore, is simply a red-shifted version 
of the spectrum they formed with, which would throw light on the physical conditions at 
the epoch of the early universe. 

Interferometric, non-resonant, detectors of gravitational radiation (with frequency con- 
tent < / < /„) use a coherent train of electromagnetic waves (of nominal frequency 
t'o ^ fu) folded into several beams, and at one or more points where these intersect, mon- 
itor relative fluctuations of frequency or phase (homodyne detection). The observed low 
frequency fluctuations are due to several causes: (a) frequency variations of the source of 
the electromagnetic signal about z/q, (b) relative motions of the electromagnetic source and 
the mirrors (or amplifying transponders) that do the folding, (c) temporal variations of the 
index of refraction along the beams, and, (d) according to general relativity, to any time- 
variable gravitational fields present, such as the transverse-traceless metric curvature of a 
passing plane gravitational wave train. To observe gravitational waves in this way, it is 
thus necessary to control, or monitor, the other sources of relative frequency fluctuations, 
and, in the data analysis, to use optimal algorithms based on the different characteristic 
interferometer responses to gravitational waves (the signal) and to the other sources (the 
noise) P]. By comparing phases of electromagnetic beams referenced to the same frequency 
generator and propagated along non-parallel equal-length arms, frequency fluctuations of 
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the frequency reference can be removed and gravitational wave signals at levels many orders 
of magnitude lower can be detected. 

In the present single-spacecraft Doppler tracking observations, for instance, many of the 
noise sources can be either reduced or calibrated by implementing appropriate microwave 
frequency links and by using specialized electronics 4], so the fundamental limitation is 
imposed by the frequency (time-keeping) fluctuations inherent to the reference clock that 
controls the microwave system. Hydrogen maser clocks, currently used in Doppler tracking 
experiments, achieve their best performance at about 1000 seconds integration time, with 
a fractional frequency stability of a few parts in 10~^^. This is the reason why these one- 
arm interferometers in space (which have one Doppler readout and a "3-pulse" response to 
gravitational waves |^) are most sensitive to millihertz gravitational waves. This integration 
time is also comparable to the microwave propagation (or "storage") time 2L/c to spacecraft 
en route to the outer solar system (for example L ^ 5 — 8 AU for the Cassini spacecraft) Qj. 

Next-generation low-frequency interferometric gravitational wave detectors in solar orbits, 
such as the Laser Interferometer Space Antenna (LISA) mission [7], have been proposed to 
achieve greater sensitivity to millihertz gravitational waves. Since the armlengths of these 
space-based interferometers can differ by few percent, the direct recombination of the two 
beams at a photo detector will not however effectively remove the laser frequency noise. This 
is because the frequency fluctuations of the laser will be delayed by different amounts within 
the two unequal length arms. In order to cancel the laser frequency noise, the time-varying 
Doppler data must be recorded and post-processed to allow for arm- length differences 
The data streams will have temporal structure, which can be described as due to many- 
pulse responses to (5-function excitations, depending on time-of-fiight delays in the response 
functions of the instrumental Doppler noises and in the response to incident plane-parallel, 
transverse, and traceless gravitational waves. 

LISA will consists of three spacecraft orbiting the sun. Each spacecraft will be equipped 
with two lasers sending beams to the other two (~0.03 AU away) while simultaneously 
measuring the beat frequencies between the local laser and the laser beams received from the 
other two spacecraft. The analysis of TDI presented in this article will assume a successful 
prior removal of any first-order Doppler beat notes due to relative motions giving six 
residual Doppler time series as the raw data of a stationary time delay space interferometer. 
Following we will regard LISA not as constituting one or more conventional 



4 



Michelson interferometers, but rather, in a symmetrical way, a closed array of six one-arm 
delay lines between the test masses. In this way, during the course of the article, we will 
show that it is possible to synthesize new data combinations that cancel laser frequency 
noises, and estimate achievable sensitivities of these combinations in terms of the separate 



and relative 
fcf. 



y simple single arm-responses both to gravitational wave and instrumental noise 



m, m)- 

In contrast to Earth-based interferometers, which operate in the long-wavelength limit 
(LWL) (arm lengths << gravitational wavelength ~c//o, where /o is a characteristic fre- 
quency of the GW), LISA will not operate in the LWL over much of its frequency band. 
When the physical scale of a free mass optical interferometer intended to detect gravitational 
waves is comparable to or larger than the GW wavelength, time delays in the response of the 
instrument to the waves, and travel times along beams in the instrument, cannot be ignored 
and must be allowed for in computing the detector response used for data interpretation. It 
is convenient to formulate the instrumental responses in terms of observed differential fre- 
quency shifts — for short, Doppler shifts — rather than in terms of phase shifts usually used 
in interferometry, although of course these data, as functions of time, are interconvertible. 

This first review article on TDI is organized as follows. In Section |n] we provide an 
overview of the physical and historical motivations of TDI. In Section UTTl we summarize the 
one-arm Doppler transfer functions of an optical beam between two carefully shielded test 
masses inside each spacecraft resulting from (i) frequency fluctuations of the lasers used in 
transmission and reception, (ii) fluctuations due to non-inertial motions of the spacecraft, 
(iii) beam-pointing fluctuations and shot noise '13]. Among these, the dominant noise is from 
the frequency fluctuations of the lasers and is several orders (perhaps 7 or 8) above the other 
noises. This noise must be very precisely removed from the data in order to achieve the GW 
sensitivity at the level set by the remaining Doppler noise sources which are at a much lower 
level and which constitute the noise floor after the laser frequency noise is suppressed. We 
show that this can be accomplished by shifting and linearly combining the twelve one-way 
Doppler data LISA will measure. The actual procedure can easily be understood in terms 
of properly defined time-delay operators that act on the one-way Doppler measurements. 
We develop a formalism involving the algebra of the time-delay operators which is based on 
the theory of rings and modules and computational commutative algebra. We show that 
the space of all possible interferometric combinations cancelling the laser frequency noise is 
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a module over the polynomial ring in which the time-delay operators play the role of the 



indeterminates. The module, in the literature, is called the Module of Syzygies [1^. We 
show that the module can be generated from four generators, so that any data combination 
cancelling the laser frequency noise is simply a linear combination formed from these gener- 
ators. We would like to emphasize that this is the mathematical structure underlying TDI 
in LISA. 

In Section IIVI specific interferometric combinations are then derived, and their physical 
interpretations are discussed. The expressions for the Sagnac interferometric combinations, 
{a, f3, 7, () are first obtained; in particular, the symmetric Sagnac combination (, for which 
each raw data set needs to be delayed by only a single arm transit time, distinguishes itself 
against all the other TDI combinations by having a higher order response to gravitational 
radiation in the LWL when the spacecraft separations are equal. We then express the 
Unequal-arm Michelson combinations, {X, Y, Z), in terms of the a, f3, 7, and ( combinations 
with further transit time delays. One of these interferometric data combinations would still 
be available if the links between one pair of spacecraft were lost. Other TDI combinations, 
which rely on only four of the possible six inter-spacecraft Doppler measurements (denoted 
P, E and U) are also presented. They would of course be quite useful in case of potential 
loss of any two inter-spacecraft Doppler measurements. 

Time-Delay Interferometry so formulated presumes the spacecraft-to-spacecraft light- 
travel-times to be constant in time, and independent from being up- or down-links. Re- 
duction of data from moving interferometric laser arrays in solar orbit will in fact encounter 
non-symmetric up- and downlink light time differences that are significant, and need to be 
accounted for in order to exactly cancel the laser frequency fluctuations 
Section 13 we show that, by introducing a set of non-commuting time-delay operators, there 
exists a quite general procedure for deriving generalized TDI combinations that account for 
the effects of time-dependence of the arms. Using this approach it is possible to derive "fiex- 
free" expression for the unequal-arm Michelson combinations Xi, and obtain the generalized 
expressions for all the TDI combinations jl?! . 



In Section IVII we address the question of maximization of the LISA signal-to-noise-ratio 
(SNR) to any gravitational wave signal present in its data. This is done by treating the SNR 
as a functional over the space of all possible TDI combinations. As a simple application of 
the general formula we have derived, we apply our results to the case of sinusoidal signals 
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randomly polarized and randomly distributed on the celestial sphere. We find that the 
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standard LISA sensitivity figure derived for a single Michelson Interferometer 
can be improved by a factor of \/2 in the low-part of the frequency band, and by more than 
in the remaining part of the accessible band. Further, we also show that if the location of 
the GW source is known, then as the source appears to move in the LISA reference frame, it 



is possible to optimally track the source, 
during the course of its trajectory [3], 



jy appropriately changing the data combinations 
|. As an example of such type of source, we 



consider known binaries within our own galaxy. 

This first version of our "Living Review" article on TDl does not include all the results 
of more practical and experimental nature, as well as all the aspects of TDI that the data 
analysts will need to account for when analyzing the LISA TDI data combinations. Our 
forthcoming "second edition" of this review paper will include these topics. It is worth 
mentioning that, as of today, the LISA project has endorsed TDI as its baseline technique 
for achieving the desired sensitivity to gravitational radiation. Several experimental verifi- 
cations and tests of TDI are being, and will be, performed at the NASA and ESA LISA 
laboratories. Although significant theoretical and experimental work has already been done 
for understanding and overcoming practical problems related to the implementation of TDI, 
more work on both sides of the Atlantic is still needed. Results of this undergoing effort will 
be included in the second edition of this living document. 



II. PHYSICAL AND HISTORICAL MOTIVATIONS OF TDI 

Equal-arm interferometer detectors of gravitational waves can observe gravitational ra- 
diation by cancelling the laser frequency fluctuations affecting the light injected into their 
arms. This is done by comparing phases of split beams propagated along the equal (but 
non-parallel) arms of the detector. The laser frequency fluctuations affecting the two beams 
experience the same delay within the two equal-length arms and cancel out at the pho- 
todetector where relative phases are measured. This way gravitational wave signals of di- 
mensionless amplitude less than 10"^*^ can be observed when using lasers whose frequency 
stability can be as large as roughly a few parts in 10~^^. 

If the arms of the interferometer have different lengths, however, the exact cancellation 
of the laser frequency fluctuations, say C{t), will no longer take place at the photodetector. 
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In fact, the larger the difference between the two arms, the larger will be the magnitude of 
the laser frequency fluctuations affecting the detector response. If Li and L2 are the lengths 
of the two arms, it is easy to see that the amount of laser relative frequency fluctuations 
remaining in the response is equal to (units in which the speed of light c = 1) 

AC{t) =C{t-2Li) -C{t-2L2) . (1) 

In the case of a space-based interferometer such as LISA, whose lasers are expected to display 
relative frequency fluctuations equal to about 10~^^/ \J Hz in the millihertz band, and whose 
arms will differ by a few percent 0, equation (PJ implies the following expression for the 
amplitude of the Fourier components of the uncanceled laser frequency fluctuations (an over 
imposed tilde denotes the operation of Fourier transform) 

|AC'(/)|^|C(/)|47r/|(Li-L2)| . (2) 

At / = 10^'^ Hz, for instance, and assuming |Li — L2I — 0.5 sec, the uncanceled fluctuations 
from the laser are equal to 6.3 x lO^^^/V^Hz. Since the LISA sensitivity goal is about 
]^g-2o^ V^Hz in this part of the frequency band, it is clear that an alternative experimental 
approach for canceling the laser frequency fluctuations is needec . 

A first attempt to solve this problem was presented in 21, 
proposed there can be understood through Figured In this idealized model the two beams 
exiting the two arms are not made to interfere at a common photodetector. Rather, each 
is made to interfere with the incoming light from the laser at a photodetector, decoupling 
in this way the phase fluctuations experienced by the two beams in the two arms. Now 
two Doppler measurements are available in digital form, and the problem now becomes one 
of identifying an algorithm for digitally cancelling the laser frequency fluctuations from a 
resulting new data combination. 



22, and the scheme 



The algorithm they first proposed, and refined subsequently in j2^, required processing 
the two Doppler measurements, say yi(t) and y2^\ in the Fourier domain. If we denote 
with hi{t), h2{t) the gravitational wave signals entering into the Doppler data y2 respec- 
tively, and with rii, n2 any other remaining noise affecting yi and y2 respectively, then the 
expressions for the Doppler observables yi, 2/2 can be written in the following form 

yr{t) = C{t - 2Li) - C{t) + h,{t) + m{t) , (3) 
y2{t) = C{t - 2L2) - C{t) + h2{t) + n2{t) . (4) 
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FIG. 1: Light from a laser is split into two beams, each injected into an arm formed by pairs of 
free-falling mirrors. Since the length of the two arms, Li and L2, are different, now the light beams 
from the two arms are not recombined at one photo detector. Instead each is separately made to 
interfere with the light that is injected into the arms. Two distinct photo detectors are now used, 
and phase (or frequency) fluctuations are then monitored and recorded there. 

From Eqs. (jHl E} it is important to note the characteristic time signature of the random 
process C(t) in the Doppler responses yi , ?/2- The time signature of the noise C(t) in yi(t), 
for instance, can be understood by observing that the frequency of the signal received at 
time t contains laser frequency fluctuations transmitted 2Li seconds earlier. By subtracting 
from the frequency of the received signal the frequency of the signal transmitted at time t, 
we also subtract the frequency fluctuations C{t) with the net result shown in Eg. (IHl) . 
The algorithm for cancelling the laser noise in the Fourier domain suggested in [2l| works 
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as follow. If we take an infinitely long Fourier transform of the data yi, the resulting 
expression of yi in the Fourier domain becomes (see eq. ©) 

mif) = C{f) [e^-^^^^ -i]+\{f) + n,{f) . (5) 

If the arm length Li is known exactly, we can use the yi data to estimate the laser frequency 
fluctuations C{f). This can be done by dividing yi by the transfer function of the laser noise 
C into the observable yi itself. By then further multiplying yi/[e'^'^'^-^^'^ — 1] by the transfer 
function of the laser noise into the other observable y2, i.e., [e^'^*-^'^^ _ g^^Ld then subtract 
the resulting expression from y2 one accomplishes the cancellation of the laser frequency 
fluctuations. 

The problem with this procedure is the underlying assumption of being able to take an 
infinitely long Fourier transform of the data. Even if one neglects the variation in time 
of the LISA arms, by taking a finite length Fourier transform of, say, yi{t) over a time 
interval T, the resulting transfer function of the laser noise C into yi no longer will be equal 
to [e'^'^^-^'^i — 1]. This can be seen by writing the expression of the finite length Fourier 
transform of yi in the following way 

/+T r+oo 
y,{t) e'-'f' dt = / y,{t) H{t) e^-^f'dt , (6) 
T J~oo 

where we have denoted with H{t) the function that is equal to 1 in the interval [— T, +T], 
and zero everywhere else. Equation © implies that the finite-length Fourier transform y^ of 
yi{t) is equal to the convolution in the Fourier domain of the infinitely long Fourier transform 



of yii with the Fourier transform of H{t) [25| (i.e. the "Sine Function" of width 1/T). 
The key point here is that we can no longer use the transfer function [e^'^*-^-^' — 1] , z = 1, 2, 
for estimating the laser noise fluctuations from one of the measured Doppler data, without 
retaining residual laser noise into the combination of the two Doppler data yi, t/2 valid in 
the case of infinite integration time. The amount of residual laser noise remaining in the 
Fourier-based combination described above, as a function of the integration time T and type 
of "window function" used, was derived in the appendix of [8]. There it was shown that, 
in order to suppress the residual laser noise below the LISA sensitivity level identified by 
secondary noises (such as proof-mass and optical path noises) with the use of the Fourier- 
based algorithm an integration time of about six months was needed. 

A solution to this problem was suggested in [sj, which works entirely in the time-domain. 
From Eqs.(I21IH) we may notice that, by taking the difference of the two Doppler data ?/i(t), 
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y2{t), the frequency fluctuations of the laser now enter into this new data set in the following 

way 

yiit) - y2{t) = C{t - 2Li) - C{t - 2L2) + Kit) - h^it) + mit) - n^it) . (7) 

If we now compare how the laser frequency fluctuations enter into Eq. ((71) against how they 
appear in Eqs. (jHl E]) we can further make the following observation. If we time-shift the 
data yi{t) by the round trip light time in arm 2, yi{t — 2L2), and subtract from it the data 
y2(t) after it has been time-shifted by the round trip light time in arm 1, y2{t — 2Li), we 
obtain the following data set 

yi{t - 2L2) - y2{t - 2Li) = C{t - 2L,) - C{t - 2L2) + h{t - 2L2) - h2{t - 2L{) 

+ni{t - 2L2) - n2{t - 2Li) . (8) 

In other words, the laser frequency fluctuations enter into yi{t) — y2{t), and yi{t — 2L2) — 
y2{t — 2Li) with the same time structure. This implies that, by subtracting Eq. (jH|) from Eq. 
([7j) we can generate a new data set that does not contain the laser frequency fluctuations 

C{t) 

X = [y,{t) - y2{t)] - Mt - 2L2) - y2{t - 2L,)] . (9) 

The expression above of the X combination shows that it is possible to cancel the laser 
frequency noise in the time domain by properly time-shifting and linearly combining Doppler 
measurements recorded by different Doppler readouts. This in essence is what Time-Delay 
Interferometry (TDI) amounts to. In the following sections we will further elaborate and 
generalize TDI to the realistic LISA conflguration. 



III. TIME-DELAY INTERFEROMETRY 



The description of TDI for LISA is greatly simplifled if we adopt the notation shown 
in Figure El , where the overall geometry of the LISA detector is deflned. There are three 
spacecraft, six optical benches, six lasers, six proof-masses and twelve photodetectors. There 
are also six phase difference data going clock-wise and counter-clockwise around the LISA 
triangle. For the moment we will make the simplifying assumption that the array is station- 
ary, i.e. the back and fort h optic al paths between pairs of spacecraft are simply equal to 
their relative distances 



fort h op tical pc 
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FIG. 2: Schematic LISA configuration. The spacecraft are labeled 1, 2, and 3. The optical paths 
are denoted by Li, L[ where the index i corresponds to the opposite spacecraft. The unit vectors 
fij point between pairs of spacecraft, with the orientation indicated. 

Several notations have been used in this context. The double index notation recently 
employed in 16] where six quantities are involved is self-evident. However, when algebraic 
manipulations are involved the following notation seems more convenient to use. The space- 
craft are labeled 1, 2, 3 and their separating distances are denoted Li, L2, L3, with Lj being 
opposite spacecraft i. We orient the vertices 1,2,3 clockwise in figure |21 Unit vectors be- 
tween spacecraft are hi, oriented as indicated in figure El We index the phase difference data 
to be analyzed as follows: The beam arriving at spacecraft i has subscript i and is primed or 
unprimed depending on whether the beam is traveling clockwise or counter-clockwise (the 
sense defined here with reference to figure 12} around the LISA triangle respectively. Thus, 
as seen from the figure, Si is the phase difference time series measured at reception at space- 
craft 1 with transmission from spacecraft 2 (along L3). Similarly, s'^ is the phase difference 
series derived from reception at spacecraft 1 with transmission from spacecraft 3. The other 
four one-way phase difference time series from signals exchanged between the spacecraft are 
obtained by cyclic permutation of the indices: 1 ^ 2 ^ 3 ^ 1. We also adopt a notation 
for delayed data streams, which will be convenient later for algebraic manipulations. We 
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define the three time-delay operators Pj, z = 1, 2, 3 where for any data stream x{t), 



Vix{t) = x{t - Li) , (10) 

where Li,i = 1,2,3 are the light travel times along the three arms of the LISA trian- 
gle (the speed of light c is assumed to be unity in this article). Thus, for example, 
V2Si{t) = si(t - La), V2V^si{t) = Slit - L2- L3) = V^V2Si{t), etc. Note that the 
operators commute here. This is because the arm-lengths have been assumed to be constant 
in time. If the Li are functions of time then the operators no longer commute 3|, 
as will be described in section IV. Six more phase difference series result from laser beams 
exchanged between adjacent optical benches within each spacecraft; these are similarly in- 
dexed as Ti, r/, 2 = 1,2, 3. The proof-mass-plus-optical-bench assemblies for LISA spacecraft 
number 1 are shown schematically in figure 2. The photo receivers that generate the data 
Si, s'l, Ti, and t[ at spacecraft 1 are shown. The phase fluctuations from the six lasers, 
which need to be canceled, can be represented by six random processes where Pi^p'^ 

are the phases of the lasers in spacecraft i on the left and right optical benches respectively 
as shown in the figure. Note that this notation is in the same spirit as in the references 
QQi„wKicK.ov,„g spacecraft a„a,s have been ana,v.ed. 

We extend the cyclic terminology so that at vertex i (i = 1, 2, 3) the random displacement 
vectors of the two proof masses are respectively denoted by 5i{t),5[{t), and the random 
displacements (perhaps several orders of magnitude greater) of their optical benches are 
correspondingly denoted by Aj(t), ^[{t) where the primed and unprimed indices correspond 
to the right and left optical benches respectively. As pointed out in 13], the analysis does 
not assume that pairs of optical benches are rigidly connected, i.e. Aj 7^ A^, in general. The 
present LISA design shows optical fibers transmitting signals both ways between adjacent 
benches. We ignore time-delay effects for these signals and will simply denote by yUi(t) the 
phase fluctuations upon transmission through the fibers of the laser beams with frequencies 
Vi, and v[. The ^i{t) phase shifts within a given spacecraft might not be the same for 
large frequency differences Vi — v[. For the envisioned frequency differences (a few hundred 
meg ahertz), however, the remaining fluctuations due to the optical fiber can be neglected 
|l3l |. It is also assumed that the phase noise added by the fibers is independent of the 
direction of light propagation through them. For ease of presentation, in what follows we 
will assume the center frequencies of the lasers to be the same, and denote this frequency 
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FIG. 3: Schematic diagram of proof-masses-plus-optical-benches for a LISA spacecraft. The left- 
hand bench reads out the phase signals si and ri. The right hand bench analogously reads out s'l 
and r{. The random displacements of the two proof masses and two optical benches are indicated 
(lower case 6i,d'j^ for the proof masses, upper case Aj, for the optical benches). 

by uq. 

The laser phase noise in Sg is therefore equal to T>iP2{t) — p'^{t). Similarly, since S2 is 
the phase shift measured on arrival at spacecraft 2 along arm 1 of a signal transmitted 
from spacecraft 3, the laser phase noises enter into it with the following time signature: 
'Dip'^{t) — P2{t)- Figure 2 endeavors to make the detailed light paths for these observations 
clear. An outgoing light beam transmitted to a distant spacecraft is routed from the laser 
on the local optical bench using mirrors and beam splitters; this beam does not interact 
with the local proof mass. Conversely, an incoming light beam from a distant spacecraft 
is bounced off the local proof mass before being reflected onto the photo receiver where it 
is mixed with light from the laser on that same optical bench. The inter-spacecraft phase 
data are denoted Si and s'-^ in figure 2. Beams between adjacent optical benches within a 
single spacecraft are bounced off proof masses in the opposite way. Light to be transmitted 
from the laser on an optical bench is first bounced off the proof mass it encloses and then 
directed to the other optical bench. Upon reception it does not interact with the proof mass 
there, but is directly mixed with local laser light, and again down converted. These data 
are denoted Ti and t[ in figure 2. 
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The expressions for the Sj,s- and rj,r/ phase measurements can now be developed from 
figures 1 and 2, and they are for the particular LISA configuration in which all the lasers have 
the same nominal frequency z/q, and the spacecraft are stationary with respect to each other. 
Consider the s[{t) process (equation (jlHj) below). The photo receiver on the right bench of 
spacecraft 1, which (in the spacecraft frame) experiences a time- varying displacement A'^^, 
measures the phase difference s'l by first mixing the beam from the distant optical bench 
3 in direction h2, and laser phase noise and optical bench motion A3 that have been 
delayed by propagation along L2, after one bounce off the proof mass {6[), with the local 
laser light (with phase noise p[). Since for this simplified configuration no frequency offsets 
are present, there is of course no need for any heterodyne conversion j^. 

In equation ()12|1 the ri measurement results from light originating at the right-bench laser 
{p[, bounced once off the right proof mass and directed through the fiber (incurring 
phase shift /ii()f:)), to the left bench, where it is mixed with laser light (pi). Similarly the 
right bench records the phase differences s'l and t[. The laser noises, the gravitational wave 
signals, the optical path noises, and proof-mass and bench noises, enter into the four data 
streams recorded at vertex 1 according to the following expressions jl^j ]: 



n 



^opt. path 



T^ZP2 -Pl + ^Q -2^3 • 5i + 723 ■ Ai + ^3 ■ V^K'^ 



p\ - pi - 2 z/Q n2 ■ {6[ - A'l) + /ii . 



gw 



'opt. path 



+ T^iVz - p'l + t'o 2^2 ■ 5[-n2- A'^ - n2 ■ 2^2 A3 



pi - pi + 2 uq hs ■ {61 - Ai) + fii 



(11) 
(12) 

(13) 

(14) 



Eight other relations, for the readouts at vertices 2 and 3, are given by cyclic permutation 
of the indices in equations (fTT|) - (fT^ . 



The gravitational wave phase signal components, sf",s^™ 



1, 2, 3, in equations 



(fTT|) and are given by integrating with respect to time the equations (1), and (2) of 
reference ll|, which relate metric perturbations to optical frequency shifts. The optical 
path phase noise contributions, s°^^' P^^^^g°P^- p^^^^ which include shot noise from the low 
signal-to-noise ratio (SNR) in the links between the distant spacecraft, can be derived from 
the corresponding term given in The Tj, r/ measurements will be made with high SNR 
so that for them the shot noise is negligible. 
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IV. ALBEGRAIC APPROACH TO CANCELLING LASER AND OPTICAL 
BENCH NOISES 



In ground based detectors the arms are chosen to be of equal length so that the laser light 
experiences identical delay in each arm of the interferometer. This arrangement precisely 
cancels the laser frequency/phase noise at the photodetector. The required sensitivity of the 
instrument can thus only be achieved by near exact cancellation of the laser frequency noise. 
However, in LISA it is impossible to achieve equal distances between spacecraft and the laser 
noise cannot be cancelled in this way. It is possible to combine the recorded data linearly 
with suitable time-delays corresponding to the three arm-lengths of the giant triangular 
interferometer so that the laser phase noise is cancelled. Here we present a systematic 
method based on modules over polynomial rings which guarantees all the data combinations 
that cancel both the laser phase and the optical bench motion noises. 

We first consider the simpler case, where we ignore the optical-bench motion noise and 
consider only the laser phase noise. We do this because the algebra is somewhat simpler 
and the method is easy to apply. The simplification amounts to physically considering each 
spacecraft rigidly carrying the assembly of lasers, beam-splitters and photodetectors. The 
two lasers on each spacecraft could be considered to be locked, so effectively there would be 
only one laser on each spacecraft. This mathematically amounts to setting Aj = = and 
Pi — Pi- The scheme we describe here for laser phase noise can be extended in a straight 
forward way to include optical bench motion noise, which we address in the last part of this 
section. 

The data combinations, when only the laser phase noise is considered, consist of the 
six suitably delayed data streams (inter-spacecraft), the delays being integer multiples of 
the light travel times between spacecraft, which can be conveniently expressed in terms of 
polynomials in the three delay operators Vi, V2, V^. The laser noise cancellation condition 
puts three constraints on the six polynomials of the delay operators corresponding to the 
six data streams. The problem therefore consists of finding six tuples of polynomials which 
satisfy the laser noise cancellation constraints. These polynomial tuples form a module ^ 

^ A module is an abelian group over a ring as contrasted with a vector space which is an abelian group over 
a field. The scalars form a ring and just like in a veetor space, scalar multiplication is defined. However, in 
a ring the multiplicative inverses do not exist in general for the elements, which makes all the difference! 



16 



called in the literature, the module of syzygies. There exist standard methods for obtaining 
the module, by that we mean, methods for obtaining the generators of the module so that 
the linear combinations of the generators generate the entire module. The procedure first 
consists of obtaining a Grobner basis for the ideal generated by the coefficients appearing 
in the constraints. This ideal is in the polynomial ring in the variables T>i^T>2-,T>^ over the 
domain of rational numbers (or integers if one gets rid of the denominators). To obtain the 
Grobner basis for the ideal, one may use the Buchberger algorithm or use an application such 



as Mathematica 



2q| . From the Grobner basis there is a standard way to obtain a generating 



set for the required module. This procedure has been described in the literature 
We thus obtain seven generators for the module. However, the method does not guarantee a 
minimal set and we find that a generating set of 4 polynomial six tuples suffice to generate 
the required module. Alternatively, we can obtain generating sets by using the software 
Macaulay 2. 

The importance of obtaining more data combinations is evident: they provide the neces- 
sary redundancy - different data combinations produce different transfer functions for GW 
and the system noises so specific data combinations could be optimal for given astrophysi- 
cal source parameters in the context of maximizing SNR, detection probability, improving 
parameter estimates etc. 



A. Cancellation of laser phase noise 

We now only have six data streams: Si and s'^ where, i = 1,2, 3. These can be regarded 
as 3 component vectors s and s' respectively. The six data streams with terms containing 
only the laser frequency noise are: 

Si = V3P2 -pi, 

s'l = T)2P3-Pu (15) 

and their cyclic permutations. 

Note that we have excluded intentionally from the data, additional phase fiuctuations 
due to the GW signal, and noises such as the optical-path noise, proof-mass noise etc. Since 
our immediate goal is to cancel the laser frequency noise we have only kept the relevant 
terms. Combining the streams for cancelling the laser frequency noise will introduce transfer 
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functions for the other noises and the GW signal. This is important and will be discussed 
subsequently in the article. 

The goal of the analysis is to add suitably delayed beams together so that the laser frequency 
noise terms add up to zero. This amounts to seeking data combinations that cancel the laser 
frequency noise. In the notation/formalism that we have invoked, the delay is obtained 
by applying the operators Vk to the beams Si and s[. A delay of kiLi + /c2-Z^2 + ksL^ is 
represented by the operator 'D^^'D^^'D^^ acting on the data, where /ci, ^2 and are integers. 
In general a polynomial in P^, which is a polynomial in three variables, applied to say Si 
combines the same data stream si{t) with different time-delays of the form kiLi + k2L2 + 
^3/73. This notation conveniently rephrases the problem. One must find six polynomials say 

g,(Pi,D2,^^3), g:(Pi,p2,2^3), z = 1, 2, 3 such that: 

3 

^ qiSi + q[s'i = . (16) 

i=l 

The zero on the R.H.S. of the above equation signifies zero laser phase noise. 

It is useful to express Eq. fjl5|) in matrix form. This allows us to obtain a matrix operator 
equation whose solutions are q and q' where and q^ are written as column vectors. We 
can similarly express Si, s'^,pi as column vectors s, s', p respectively. In matrix form Eq. (|TH|l 
become: 



s 



■ p , s' = D ■ p , (17) 



where, D is a 3 x 3 matrix given by, 

^ -1 I?2 ^ 

D= P3 -1 . (18) 

The exponent 'T' represents the transpose of the matrix. Eq. becomes: 

q^ • s + q'^ ■ s' = (q^ ■ + q'^ ■ D) ■ p = . (19) 

where we have taken care to put p on the right-hand-side of the operators. Since the above 
equation must be satisfied for an arbitrary vector p, we obtain a matrix equation for the 
polynomials (q, q'): 

q^ ■ + q' ■ D = . (20) 

Note that since the Vk commute, the order in writing these operators is unimportant. In 
mathematical terms, the polynomials form a commutative ring. 
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B. Cancellation of laser phase noise in the unequal arm interferometer 



The use of commutative algebra is very conveniently illustrated with the help of the sim- 
pler example of the unequal arm interferometer. Here there are only two arms instead of 
three as we have for LISA and the mathematics is much simpler and so it easy to see both 
physically and mathematically how commutative algebra can be applied to this problem of 
laser phase noise cancellation. The procedure is well known for the unequal arm interferom- 
eter, but here we will describe the same method but in terms of the delay opertors that we 
have introduced. 

Let (j){t) denote the laser phase noise entering the laser cavity as shown in Fig. 0J Consider 
this light (f){t) making a round trip around arm 1 whose length we take to be Li. If we 
interfere this phase with the incoming light we get the phase (f)i{t), where, 

0i(t) = 0(t - 2Li) - = {Vl - (21) 

The second expression we have written in terms of the delay operators. This makes the 
procedure transparent as we shall see. We can do the same for the arm 2 to get another 
phase 4>2(t), where, 

Mt) = <Pit - 2L2) - 0(t) = {Vl - (22) 

Clearly, if Li 7^ L2, then the difference in phase 02 (^) ~ 0i(^) is not zero and the laser phase 
noise does not cancel out. However, if one further delays the phases 0i(t) and 02 (^) and 
constructs the following combination: 

X{t) = [02(t - 2Li) - 02(t)] - [01 (t - 2L2) - 0i(t)], (23) 

then the laser phase noise does cancel out. We have already encountered this combination at 
;he end of section II, Eq.ijH)), when it was first proposed by Tinto and Armstrong in reference 
^. The cancellation of laser frequency noise becomes obvious from the operator algebra in 

the following way. In the operator notation: 

x{t) = (vl-i)Mt)-i'Dl-i)Mt) 

= [{Vj - l){Vl - 1) - {Vl - l){Vl - l)]0(t) 

= 0. (24) 

From this one immediately sees that just the commutativity of the operators has been used 
to cancel the laser phase noise. The basic idea was to compute the lowest common multiple 
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FIG. 4: Schematic diagram of the unequal-arm Michelson interferometer. The beam shown cor- 
responds to the term (Vf — 1){T>2 — l)'^(i) in X{t) which is first sent around arm 1 fohowed by 
arm 2. The second beam (not shown) is first sent around arm 2 and then through arm 1. The 
difference in these two beams constitutes X{t). 

(L.C.M.) of the polynomials Vf — l and P| — 1 (in this case the L.C.M. is just the product, 
because the polynomials are relatively prime) and use this fact to construct X{t) in which 
the laser phase noise is cancelled. The operation is shown physically in Fig. 0] The notions 
of commutativity of polynomials, L.C.M. etc. belong to the field of commutative algebra. In 
fact we will be using the notion of a Grobner basis which is in a sense the generalization of 
the notion of the greatest common divisor or in short gcd. Since LISA has three spacecraft 
and six inter-spacecraft beams, the problem of the unequal arm interferometer only gets 
technically more complex; in principle the problem is the same as in this simpler case. Thus 
the simple operations which were performed here to obtain a laser noise free combination 
X{t) are not sufficient and more sophisticated methods need to be adopted from the field 
of commutative algebra. We address this problem in the forthcoming text. 
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C. The Module of Syzygies 



Eq. ()20|) has non-trivial solutions. Several solutions have been exhibited in ll|,ll3|. We 
merely mention these solutions here; in the next section we will discuss them in detail. 
The solution ( is given by — q-^ = q'^ = (Pi, 1^2, ^^s)- The solution a is described by 
q-^ = —(1,1)3,1)11)3) and q'"'" = (1, D1D2J ^^2)- The solutions P and 7 are obtained from 
a by cyclically permuting the indices of T'fc,q and q'. These solutions are important, be- 
cause they consist of polynomials with lowest possible degrees and thus are simple. Other 
solutions containing higher degree polynomials can be generated conveniently from these 
solutions. Since the system of equations is linear, linear combinations of these solutions are 
also solutions to Eq. (^0]) . 

However, it is important to realize that we do not have a vector space here. Three 
independent constraints on a six tuple do not produce a space which is necessarily generated 
by three basis elements. This conclusion would follow if the solutions formed a vector 
space but they do not. The polynomial six-tuple q, q' can be multiplied by polynomials in 
T'i,T'2)'P3 (scalars) which do not form a field. So that the inverse in general does not exist 
within the ring of polynomials. We therefore have a module over the ring of polynomials in 
the three variables Di,D2,^^3- First we present the general methodology for obtaining the 
solutions to (pUj) and then apply it to equations (j^ . 

There are three linear constraints on the polynomials given by the equations ()20|). Since 
the equations are linear the solutions space is a submodule of the module of six-tuples of 
polynomials. The module of six-tuples is a free module, i.e. it has six basis elements that 
not only generate the module but are linearly independent. A natural choice of the basis 
is fm = (0, 1, 0) with 1 in the m-th place and everywhere else; m runs from 1 to 6. 
The definitions of generation (spanning) and linear independence are the same as that for 
vector spaces. A free module is essentially like a vector space. But our interest lies in its 
submodule which need not be free and need not have just three generators as it would seem 
if we were dealing with vector spaces. 

The problem at hand is of finding the generators of this submodule i.e. any element of 
the submodule should be expressible as a linear combination of the generating set. In this 
way the generators are capable of spanning the full submodule or generating the submodule. 
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In order to achieve our goal, we rewrite the Eq. ()20j) exphcitly component-wise: 

qi + q'i- Vsq'^ - P2g3 = , 
?2 + ^2 - ^1^3 - ^3^1 = , 

Qs + q's- '^2q[ - 'Diq2 = . (25) 
The first step is to use Gaussian ehmination to obtain qi and q2 in terms of q3,q[,q2,q3- 



qi = -q'l + ^^3^2 + ^^2^3 , 

q2 = -q'2 + T^iq'i + T^zqi 

= -Vsq[ - (1 - Vl)q'^ + VW3 + T^2V3qz , (26) 

and then substitute these values in the third equation to obtain a linear implicit relation 
between q3,qi,q2,q3- We then have: 

(1 - V{D2V,)q, + (D1D3 - V2)q[ + Pi(l - Vl)q'2 + (1 - VDq'.^ = . (27) 

Obtaining solutions to Eq. (jTTj) amounts to solving the problem since the the remaining 
polynomials gi, q2 have been expressed in terms of q3,q'i,q2, qs in fOUj) . Note that we cannot 
carry on the Gaussian elimination process any further, because none of the polynomial 
coefficients appearing in Eq.()27p have an inverse in the ring. 

We will assume that the polynomials have rational coefficients i. e the coefficients belong 
to Q the field of the rational numbers. The set of polynomials form a ring - the polynomial 
ring in three variables which we denote by 7^ = D25 ^^3]- The polynomial vector 

(?3) Q'l; 5'25 ^3) ^ The set of solutions to (jTTj) is just the kernel of the homomorphism 

(f : 7?.^ IZ where the polynomial vector (i?3, Q'l, q'2, ^3) is mapped to the polynomial 
{l-ViV2V^)q^ + {ViV^-V2)q[+Vi{l-V'j^q'^ + {l-Vl)q'^. Thus the solution space kenp 
is a submodule of 7^^. It is called the module of syzygies in the literature. The generators of 
this module can be obtained from standard methods available in the literature. We briefly 
outline the method given in the books by Becker et al. 
below. The details have been included in appendix A. 



27[ and Kreuzer and Robbiano 
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D. Grobner Basis 



The first step is to obtain the Grobner basis for the ideal U generated by the coefficients 
in Eq. (fTTjl : 

ui = l-ViV2V^, U2=ViVs-V2, U3 = V^{l-Vl), Ui = l-Vj. (28) 

The ideal 14 consists of linear combinations of the form 'Y^ViUi where fj, i = 1,...,4 are 
polynomials in the ring TZ. There can be several sets of generators for U. A Grobner basis 
is a set of generators which is 'small' in a specific sense. 

There are several ways to look at the theory of Grobner basis. One way is, suppose we are 
given polynomials gi, g2, Qm one variable over say Q and we would like to know whether 
another polynomial / belongs to the ideal generated by the g^s. A good way to decide the 
issue would be to first compute the gcd (greatest common divisor) g of gi, g2, gm and 
check whether / is a multiple of g. One can achieve this by doing the long division of / by 5^ 
and checking whether the remainder is zero. All this is possible because Q[x] is a Euclidean 
domain and also a principle ideal domain (PID) wherein any ideal is generated by a single 
element. Therefore we have essentially just one polynomial - the gcd - which generates the 
ideal generated by gi-, g2i ■ ■ ■ 1 Qm- The ring of integers or the ring of polynomials in one 
variable over any field are examples of PIDs whose ideals are generated by single elements. 
However, when we consider more general rings (not PIDs) like the one we are dealing with 
here, we do not have a single gcd but a set of several polynomials which generates an ideal 
in general. A Grobner basis of an ideal can be thought of as a generalization of the gcd. In 
the univariate case, the Grobner basis reduces to the gcd. 

Grobner basis theory generalizes these ideas to multivariate polynomials which are nei- 
ther Euclidean rings nor PIDs. Since there is in general not a single generator for an ideal, 
Grobner basis theory comes up with the idea of dividing a polynomial with a set of polyno- 
mials, the set of generators of the ideal, so that by successive divisions by the polynomials 
in this generating set of the given polynomial, the remainder becomes zero. Clearly, every 
generating set of polynomials need not possess this property. Those special generating sets 
that do possess this property (and they exist!) are called Grobner bases. In order for a divi- 
sion to be carried out in a sensible manner, an order must be put on the ring of polynomials, 
so that the final remainder after every division is strictly smaller than each of the divisors in 
the generating set. A natural order exists on the ring of integers or on the polynomial ring 
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Q{x); the degree of the polynomial decides the order in Q{x). However, even for polynomials 
in two variables there is no natural order apriori (Is x"^ + y greater or smaller than x + y^?). 
But one can, by hand as it were, put an order on such a ring by saying x >> y, where >> 
is an order, called the lexicographical order. We follow this type of order, Vi » V2 >> P3 
and ordering polynomials by considering their highest degree terms. It is possible to put 
different orderings on a given ring which then produce different Grobner bases. Clearly, 
a Grobner basis must have 'small' elements so that division is possible and every element 
of the ideal when divided by the Grobner basis elements leaves zero remainder, i.e. every 
element modulo the Grobner basis reduces to zero. 

In the literature, there exists a well-known algorithm called the the Buchberger algorithm 
which may be used to obtain the Grobner basis for a given set of polynomials in the ring. 
So a Grobner basis of U can be obtained from the generators Ui given in Eq. (PHjl using this 
algorithm. It is essentially again a generalization of the usual long division that we perform 
on univariate polynomials. More conveniently, we prefer to use the well known application 
'Mathematica'. Mathematica yields a 3 element Grobner basis Q for U: 

G = {Vl - 1, Vl - 1, Pi - . (29) 

One can easily check that all the Ui of Eq. (|28|) are linear combinations of the polynomials 
in Q and hence Q generates U. One also observes that the elements look 'small' in the 
order mentioned above. However, one can satisfy oneself that ^ is a Grobner basis by using 
the standard methods available in the literature. One method consists of computing the 
S-polynomials (see Appendix A) for all the pairs of the Grobner basis elements and checking 
whether these reduce to zero modulo Q. 

This Grobner basis of the ideal U is then used to obtain the generators for the module 
of syzygies. Note that although the Grobner basis depends on the order we choose among 
the Pfc, the module itself is independent of the order 



E. Generating Set for the Module of Syzygies 



The generating set for the module is obtained by further following the procedure in the 



literature 
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28[. The details are given in Appendix A, specifically for our case. We obtain 



7 generators for the module. These generators do not form a minimal set and there are 
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relations between them; in fact this method does not guarantee a minimum set of generators. 
These generators can be expressed as hnear combinations of 7,^ and also in terms of 
given below in Eq. (j3(jp . The importance in obtaining the 7 generators 
is that the standard theorems guarantee that these 7 generators do in fact generate the 
required module. Therefore, from this proven set of generators we can check whether a 
particular set is in fact a generating set. We present several generating sets below. 

Alternatively, we may use a software package called Macaulay 2 which directly calculates 
the generators given the the equations (|^. Using Macaulay 2, we obtain six generators. 
Again, Macaulay's algorithm does not yield a minimal set; we can express the last two 
generators in terms of the first four. Below we list this smaller set of four generators in the 
order X = (gi, gs, gs, q'v 4-1^ q's)- 

= {V2 - 0, 1 - Vl 0, V^Vs - Pi, Vl - 1) , 

= {-V,V2, -1, -Vl, Vs, 1, V2V^) . (30) 

Note that the last three generators are just X^"^^ = (, X^^^ = a, X^^^ = (3. An extra generator 
X*^^) is needed to generate all the solutions. 

Another set of generators which may be useful for further work is a Grobner basis of a 
module. The concept of a Grobner basis of an ideal can be extended to that of a Grobner 
basis of a submodule of {K[xi,X2, .■■,Xj^)^ where is a field, since a module over the 
polynomial ring can be considered as generalization of an ideal in a polynomial ring. Just 
as in the case of an ideal, a Grobner basis for a module is a generating set with spe- 
cial properties. For the module under consideration we obtain a Grobner basis using 
Macaulay 2 : 

G-W = (-Pi,-P2,-P3,Pi,P2,I?3) , 

G^^) = (V2 - V{D;, 0, 1 - vl 0, V^V^ - Pi, Vl - 1) , 

= (-1,-1)3, -PiP3,l,PiP2,I?2) , 

G'(^) = {V,{l-Vl),Vl-l,Q,Q,l-VlV,{Vl-l)) . (31) 
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Note that in this Grobner basis G^^) = C = X(2),G(2) = _ ^ _ X^^\G^^'^ = a = 

X^'^\ Only G^^^ is the new generator. 

Another set of generators are just and (. This can be checked using Macaulay 2 

or one can relate a, /5, 7 and C, to the generators X*^'^^ A= 1, 2, 3, 4 by polynomial matrices. 
In Appendix B, we express the 7 generators we obtained following the literature, in terms 
of 7 and C,. Also we express a,/3,'y and ( in terms of X^'^\ This proves that all these 
sets generate the required module of syzygies. 

The question now arises as to which set of generators we should choose which facilitates 
further analysis. The analysis is simplified if we choose a smaller number of generators. 
Also we would prefer low degree polynomials to appear in the generators so as to avoid 
cancellation of leading terms in the polynomials. By these two criteria we may choose, X^^^^ 
or a, P, 7, (. However, a, /3, 7, ( possess the additional property that this set is left invariant 
under a cyclic permutation of indices 1,2,3. It is found that this set is more convenient to 
use because of this symmetry. 

F. Canceling optical bench motion noise 

There are now twelve Doppler data streams which have to be combined in an appropriate 
manner in order to cancel the noise from the laser as well as from the motion of the optical 
benches. As in the previous case of cancelling laser phase noise, here too, we keep the relevant 
terms only, namely, those terms containing laser phase noise and optical bench motion noise. 
We then have the following expressions for the four data streams on spacecraft 1: 

si = Vs[p'2 + z/ofis ■ A2] - [pi - z/ofis ■ Ai] , (32) 

s[ = T'2[P3 - ^^0112 ■ A3] - [p[ + z/ofia ■ A[] , (33) 

n = P'l -Pi + 2uoh2 ■ A[ + /ii , (34) 

n = Pi- p'l - St'ona ■ Ai + /ii . (35) 

The other eight data streams on spacecraft 2 and 3 are obtained by cyclic permutations of 
the indices in the above equations. In order to simplify the derivation of the expressions 
cancelling the optical bench noises, we note that by subtracting eq. (j35|) from eq. (j34|) . we 
can rewriting the resulting expression (and those obtained from it by permutation of the 
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spacecraft indices) in the following form, 



zi = ]^{ti- t[) = 01 - , (36) 



where (f)[, 0i are defined as, 



(f)[ = p[ + uon2 ■ , 

01 = Pi - z^ofis • Ai , (37) 

The importance in defining these combinations is that the expressions for the data streams 
Si, s'i simplify into the following form, 

si = 1)302 - 01 ) 

s[ = ©203 - 0'i • (38) 
If we now combine the Si, s'^, and Zi in the following way, 

rii = Si- V3Z2 = ©302 - 01 , Vl' = Sl' + Zi= ©203 - 01 , (39) 
?72 = S2 - ViZ^ = Pi03 - 02 , r]2' = S2' + Z2 = 1^301 - 02 , (40) 
?73 = S3 - V2Z1 = ©201 - 03 , V3' = S3' + Z3 = Pi 02 - 03 , (41) 

we have just reduced the problem of cancelling of six laser and six optical bench noises to 
the equivalent problem of removing the three random processes, 0i, 02, and 03, from the 
six linear combinations, rji, rj^, of the one-way measurements s,, s\, and Zi. By comparing 
the equations above to equation p5|) for the simpler configuration with only three lasers, 
analyzed in the previous section, we see that they are identical in form. 



G. Physical Interpretation of the TDI combinations 

It is important to notice that the four interferometric combinations (a,/3, 7,^), which 
can be used as a basis for generating the entire TDI space, are actually synthesized Sagnac 
interferometers. This can be seen by rewriting the expression for a, for instance, in the 
following form, 

a = [r?r + V2V3' + ViV2'V2'] - [Vi + ^^3^2 + ^^12^3^2] , (42) 
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and noticing that the first square bracket on the right-hand side of equation ()42|) contains a 
combination of one-way measurements describing a hght beam propagating clockwise around 
the array, while the other terms in the second square-bracket give the equivalent of another 
beam propagating counter-clockwise around the constellation. 

Contrary to a, /?, and 7, ( can not be visualized as the difference(or interference) of two 
synthesized beams. However, it should still be regarded as a Sagnac combination since there 



exists a time-delay relationship between it and a, (3, and 7 



C - ViViVsC = Via - V2V^a + Pa" - V.J)iP + ©37 - V1V2I . (43) 

As a consequence of the time-structure of this relationship, C, has been called the Symmetrized 
Sagnac combination. 

By using the four generators, it is possible to construct several other interferometric 
combinations, such as the unequal-arm Michelson {X,Y,Z), the Beacons {P,Q,R), the 
Monitors {E,F,G), and the Relays {U,V,W). Contrary to the Sagnac combinations, these 
only use four of the six data combinations rji, rj^. For this reason they have obvious utility 
in the event of selected subsystem failures [l3|. 

These observables can be written in terms of the Sagnac observables (a, /3, 7, Q in the 
following way, 

ViX = V2V^a - V2P - Ds-f + C , 
P = C-Via, 
E = Di-PiC, 

U = V,-f-l3, (44) 



as it is easy to verify by substituting the expressions for the Sagnac combinations into the 
above equations. Their physical interpretations are schematically shown in Figure El In the 
case of the combination X, in particular, by writing it in the following form 



'igure 

B. 



X = [{r][ + V2'r]3) + V2'V2{vi + V^r]'^)] - [ivi + V3V2) + V^Vyir]', + Ps'r/s)] , (45) 
one can notice (as pointed out in and Q) 

that this combination can be visualized as the 
difference of two sums of phase measurements, each corresponding to a specific light path 
from a laser onboard spacecraft 1 having phase noise 0i. The first square-bracket term in 
equation ()45|) represents a synthesized light-beam transmitted from spacecraft 1 and made 
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FIG. 5: Schematic diagrams of the unequal-arm Michelson, Monitor, Beacon, and Relay combi- 
nations. These TDI combinations rely only on four of the six one-way Doppler measurements, as 
illustrated here. 

to bounce once at spacecraft 2 and 3 respectively. The second square-bracket term instead 
corresponds to another beam also originating from the same laser, experiencing the same 
overall delay as the first beam, but bouncing off spacecraft 3 first and then spacecraft 2. 
When they are recombined they will cancel the laser phase fluctuations exactly, having both 
experienced the same total delay (assuming stationary spacecraft). The X combinations 
should therefore be regarded as the response of a zero-area Sagnac interferometer. 



V. TIME-DELAY INTERFEROMETRY WITH MOVING SPACECRAFT 

The rotational motion of the LISA array results in a difference of the light travel times 
in the two directions around a Sagnac circuit Two time delays along each arm 

must be used, say L[ and Li for clockwise or counterclockwise propagation as they enter 
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in any of the TDI combinations. Furthermore, since and L\ not only differ from one 
another but can be time dependent (they "flex"), it was shown that the "first generation" 
TDI combinations do not completely cancel the laser phase noise (at least with present laser 
stability requirements), which can enter at a level above the secondary noises. For LISA, 
and assuming Lj ~ lOm/sec |30i|, the estimated magnitude of the remaining frequency 
fluctuations from the laser can be about 30 times larger than the level set by the secondary 
noise sources in the center of the frequency band. In order to solve this potential problem, 
it has been shown that there exist new TDI combinations that are immune to first order 
shearing (flexing, or constant rate of change of delay times). These combinations can be 
derived by using the time-delay operators formalism introduced in the previous section. 



although one has to keep in mind that now these operators no longer commute [17 1. 

In order to derive the new, "flex-free" TDI combinations we will start by taking specific 
combinations of the one-way data entering in each of the expressions derived in the previous 
section. These combinations are chosen in such a way so as to retain only one of the three 
noises 0j, z = 1,2,3 if possible. In this way we can then implement an iterative procedure 
based on the use of these basic combinations and of time-delay operators, to cancel the laser 
noises after dropping terms that are quadratic in L/c or linear in the accelerations. This 
iterative time-delay method, to first order in the velocity, is illustrated abstractly as follows. 
Given a function of time \1/ = '^{t), time delay by is now denoted either with the standard 
comma notation or by applying the delay operator Vi introduced in the previous section 

V,^ = ^, = ^{t-L,{t)) . (46) 
We then impose a second time delay Lj{t): 

VjV,^ = ^.i,- = ^{t- Lj{t) - Li{t - Lj{t))) 

- ^it-Lj{t)-Li{t) + Li{t)Lj) 

- '^,^J + ^,^JL^L, . (47) 

A third time delay Lk{t) gives: 

VkVjV,^ = ^.ijk = ^{t- Lk{t) - L^{t - Lk{t)) -L,{t- Lk{t) - Lj{t - Mt)))) 

~ + ^^ijk[Li{Lj + Lk) + LjLk] , (48) 



30 



and so on, recursively; each delay generates a first-order correction proportional to its rate 
of change times the sum of all delays coming after it in the subscripts. Commas have now 
been replaced with semicolons (l^, to remind us that we consider moving arrays. When the 
sum of these corrections to the terms of a data combination vanishes, the combination is 
called flex-free. 

Also, note that each delay operator, Vi, has a unique inverse, D~^, whose expression can 
be derived by requiring that D^^Vi = I, and neglecting quadratic and higher order velocity 
terms. Its action on a time series is 

D7'^{t) = ^{t + L,{t + L,)) . (49) 

Note that this is not like an advance operator one might expect, since it advances not by 
Li(t) but rather Li(t + Li). 



A. The Unequal- Arm Michelson 

The unequal-arm Michelson combination relies on the four measurements r^i, r^i/, 772', 
and 773. Note that the two combinations r/i + 772', 3, + ^73,2' represent the two synthesized 
two-way data measured onboard spacecraft 1, and can be written in the following form 

Til + ^2',3 = {V-iV-y -I) 01 , (50) 
Vl' + V3,2' = iT)2''D2 - I) 01 , (51) 

where / is the identity operator. Since in the stationary case any pairs of these operators 
commute, i.e. ViVji — VjiVi = 0, from equations EH) it is easy to derive the following 
expression for the unequal- arm interfero metric combination, X, which eliminates, 0i 

X = [V2'V2 - I] iVl + r/2',3) - [(^^3^^3' - /)] iVl' + V3,2'). (52) 

If, on the other hand, the time-delays depend on time, the expression of the unequal-arm 
Michelson combination above no longer cancels 0i. In order to derive the new expression 
for the unequal-arm interferometer that accounts for "flexing", let us first consider the 
following two combinations of the one-way measurements entering into the X observable 
given in equation (j^ : 

liVl' + V3;2') + iVl + ^2;3);22'] = [D2'^2^3^3' - /] 01 , (53) 
[(^1 + ^2';3) + (^1' + ^?3;2');3'3] = [£'3^3'^2'^2 " /] 01 • (54) 
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Using equations (fSHl Elj) we can use the delay technique again to finally derive the following 
expression for the new unequal-arm Michelson combination Xi that accounts for the fiexing 
effect, 

Xi = [D2D2'DyD3 - I] [{r]21 + r]i2-y) + (?731 + '7l3;2);33'] 

- [DyD3D2D2' - I] [(r^si + '7l3;2) + (^?21 + ^?12;3' ) ;2'2] • (55) 

As usual, X2 and X3 are obtained by cyclic permutation of the spacecraft indices. This ex- 
pression is readily shown to be laser-noise-free to first order of spacecraft separation velocities 
Li-, it is "fiex-free". 

B. The Sagnac Combinations 

In the above subsection we have used the same symbol X for the unequal-arm Michelson 
combination for both the rotating (i.e. constant delay times) and stationary cases. This 
emphasizes that, for this TDI combination (and, as we will see below, also for all the 
combinations including only four links) the forms of the equations do not change going from 
systems at rest to the rotating case. One needs only distinguish between the time-of-fiight 
variations in the clockwise and counter-clockwise senses (primed and unprimed delays). 

In the case of the Sagnac variables, {a, f3, 7, (), however, this is not the case as it is easy to 
understand on simple physical grounds. In the case of a for instance, light originating from 
spacecraft 1 is simultaneously sent around the array on clockwise and counterclockwise loops, 
and the two returning beams are then recombined. If the array is rotating, the two beams 
experience a different delay (the Sagnac effect), preventing the noise (pi from cancelling in 
the a combination. 

In order to find the solution to this problem let us first rewrite a in such a way to explicitly 
emphasize what it does: attempts to remove the same fiuctuations affecting two beams that 
have been made to propagated clockwise and counter-clockwise around the array, 

a = [r]r + V2'r]y + VrV2'r]2'] - [t]i + ©3^2 + ^^i^^3^2] , (56) 

where we have accounted for clockwise and counterclockwise light delays. It is straightfor- 
ward to verify that this combination no longer cancels the laser and optical bench noises. 
If, however, we expand the two terms inside the square-brackets on the right-hand-side of 
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equation (fSUj) we find that they are equal to: 

iVl' + 1^2' + 'DyV2'r]2'] = 



[Vi + 1^3V2 + V{D^iq2] = [V^V{D2 - • (57) 

If we now apply our iterative scheme to the combinations given in equation ()57|) we finally get 
the expression for the Sagnac combination, ai, that is unaffected by laser noise in presence 
of rotation, 

ai = [V^V{D2 - I] [i^v + V^'Vv + VyV2'm>] - [V^'VyVy - I][r)i + V^r^^ + V^D^r^^] . (58) 

If the delay-times are also time-dependent, we find that the residual laser noise remaining 
into the combination ai is actually equal to 

01,1231'2'3' [{U +L2 + L3) {L[ + 4 + 4) - (4 + 4 + ^3) (^1 +L2 + L3)] . (59) 

Fortunately, although first order in the relative velocities, the residual is small, as it involves 
the difference of the clockwise and counterclockwise rates of change of the propagation delays 
on the same circuit. For LISA, the remaining laser phase noises in ai, i = 1,2, 3, are several 
orders of magnitude below the secondary noises. 

In the case of (, however, the rotation of the array breaks the symmetry and therefore 
its uniqueness. However, there still exist three generalized TDI laser-noise-free data combi- 
nations that have properties very similar to (, and which can be used for the same scientific 



purposes |31||. These combinations, which we call (Ci;C25C3)) can be derived by applying 
again our time-delay operator approach. 

Let us consider the following combination of the rji, rjii measurements, each being delayed 



only once 



^73,3 - m',3 + VIA' = [D3D2 - Dy] (pi , (60) 

r]V,l - r/2,2' + '72',2' = [^3' ^2' - ^1] 01 , (61) 

where we have used the commutativity property of the delay operators in order to cancel 
the 02 and 03 terms. Since both sides of the two equations above contain only the 0i noise, 
Ci is found by the following expression: 

Ci = [D^,D2' - Di] (?73i_i/ - 7732,2 + r/12,2) - [D2D3 - Di,] (7713,3/ - 7/23,3' + ^21,1) ■ (62) 
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If the light-times in the arms are equal in the clockwise and counterclockwise senses (e.g. 
no rotation) there is no distinction between primed and unprimed delay times. In this case, 
Ci is related to our original symmetric Sagnac C by Ci = C23 — Ci- Thus for the practical 
LISA case (arm length difference < 1%), the SNR of will be the same as the SNR of 

If the delay-times also change with time, the perfect cancellation of the laser noises is no 
longer achieved in the (Ci, C2; Cs) combinations. However, it has been shown in jl^] that the 
magnitude of the residual laser noises in these combinations are significantly smaller than 
the LISA secondary system noises, making their effects entirely negligible. 

The expressions for the Monitor, Beacon, and Relay combinations, accounting for the 
rotation and flexing of the LISA array, have been derived in the literature by applying the 
time-delay iterative procedure highlighted in this section. The interested reader is referred 
to that paper for details. 

A mathematical formulation of the "second generation" TDI, which generalizes the one 
presented in Section |l3 for the stationary LISA, still needs to be derived. In the case 
when only the Sagnac effect is considered (and the delay-times remain constant in time) the 
mathematical formulation of Section IIVI can be extended in a straight forward way where 
now the six time-delays "Dj and "D- must be taken into account. The polynomial ring is 
now in these six variables and the corresponding module of syzygies can be constructed over 
this enlarged polynomial ring 3^. However, when the arms are allowed to flex, that is, the 
operators themselves are functions of time, the operators no longer commute. One must then 
resort to non-commutative Grobner basis methods. We will investigate this mathematical 
problem in the near future. 

VI. OPTIMAL LISA SENSITIVITY 



'erometric combinations have been shown to individually have rather 
isl ]. as a consequence of their different responses to gravitational ra- 



All the above inter 
different sensitivities 

diation and system noises. Since LISA has the capability of simultaneously observing a 
gravitational wave signal with many different interferometric combinations (all having dif- 
ferent antenna patterns and noises), we should no longer regard LISA as a single detector 
system but rather as an array of gravitational wave detectors working in coincidence. This 
suggests that the presently adopted LISA sensitivity could be improved by optimally com- 
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billing elements of the TDI space. 

Before proceeding with this idea, however, let us consider again the so called "second 
generation" TDI Sagnac observables: (ai, as). The expressions of the gravitational wave 
signal and the secondary noise sources entering into ai will in general be different from 



espo] 



those entering into a, the corresponding Sagnac observable derived under the assumption of 

However, the other remaining, secondary noises in LISA 



a stationary LISA array 
are so much smaller, and the rotation and systematic velocities in LISA are so intrinsically 



small, that index permutation may still be done for them [17j. It is therefore easy to derive 



the following relationship between the signal and secondary noises in ai, and those entering 



into the stationary TDI combination a 



0, 



ai(t)^a(t)-a(t-Li-L2-L3) , (63) 

where Lj , i = 1, 2, 3 are the unequal-arm lengths of the stationary LISA array. Equation 
fl63|) implies that any data analysis procedure and algorithm that will be implemented for 
the second-generation TDI combinations can actually be derived by considering the corre- 
sponding "first generation" TDI combinations. For this reason, from now on we will focus 
our attention on the gravitational wave responses of the first-generation TDI observables 
(a,/?,7,C)- 

As a consequence of these considerations, we can still regard (a, /3, 7, Q as the generators 
of the TDI space, and write the most general expression for an element of the TDI space, 
//(/), as a linear combination of the Fourier transforms of the four generators (5, /3, 7, C) 

r^(/) ^ ai(/. A) 5(/) + a2(/. A) ^(/) + 03(7, A) 7(/) + a,{f ,X) ({f) , (64) 

where the {ai{f, ^)}j=i are arbitrary complex functions of the Fourier frequency /, and of a 
vector A containing parameters characterizing the gravitational wave signal (source location 
in the sky, waveform parameters, etc.) and the noises affecting the four responses (noise 
levels, their correlations, etc.). For a given choice of the four functions {ai}f^i, rj gives an 
element of the functional space of interferometric combinations generated by {a, (3, '-/,() ■ 
Our goal is therefore to identify, for a given gravitational wave signal, the four functions 
{'^i}i=i that maximize the signal-to-noise ratio, SNR^, of the combination i], 

SNRJ = + + + df . (65) 

Jfl (hi + 02 /3n + 03 7n + a4Cnp) 
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In equation (jU^ the subscripts s and n refer to the signal and the noise parts of (5,/3,7, C) 
respectively, the angle brackets represent noise ensemble averages, and the interval of inte- 
gration (//, fu) corresponds to the frequency band accessible by LISA. 

Before proceeding with the maximization of the SNR^ we may notice from equation (jlH)) 
that the Fourier transform of the totally symmetric Sagnac combination, (, multiplied by 
the transfer function 1 — e^'^^-^^^i+'^^+^s) qq^^i be written as a linear combination of the Fourier 
transforms of the remaining three generators (5,/?, 7). Since the signal-to-noise ratio of rj 
and (1 — e^'^*'^*^-^^"'"^2+-^^))?7 are equal, we may conclude that the optimization of the signal- 
to-noise ratio of rj can be performed only on the three observables a, f3, 7. This implies the 
following redefined expression for SNR^: 



hi ttn + 02 /?n + as 7„|2) 

The SNR^ can be regarded as a functional over the space of the three complex functions 
{cii}f^i, and the particular set of complex functions that extremize it can of course be derived 
by solving the associated set of Euler-Lagrange equations. 

In order to make the derivation of the optimal SNR easier, let us first denote by x*^^^ 
and x(") the two vectors of the signals (asyPsyls) and the noises (an, Pn,'ln) respectively. 
Let us also define a to be the vector of the three functions {ai}f=i, and denote with C the 
hermitian, non-singular, correlation matrix of the vector random process x„, 

(C)., = (x(")xl")*) . (67) 

If we finally define (A)jj to be the components of the hermitian matrix x-'^-'x^*-'*, we can 
rewrite SNR^ in the following form, 

SNR? = if , (68) 

where we have adopted the usual convention of summation over repeated indices. Since the 
noise correlation matrix C is non-singular, and the integrand is positive definite or null, 
the stationary values of the signal-to-noise ratio will be attained at the stationary values of 
the integrand, which are given by solving the following set of equations (and their complex 
conjugated expressions): 

d 



da.k 



= , fc=l,2,3. (69) 
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After taking the partial derivatives, equation ()69p can be rewritten in the following form 



(C-^)..(A),,(a*), 



1,2,3 (70) 



which tells us that the stationary values of the signal-to-noise ratio of rj are equal to the 
eigenvalues of the the matrix C""*" ■ A. The result in equation (l69ll is well known in the 
theory of quadratic forms, and it is called the Rayleigh's principlej^, [s^. 

In order now to identify the eigenvalues of the matrix C"'*' ■ A, we first notice that the 
3x3 matrix A has rank 1. This implies that the matrix C^^ • A has also rank 1, as it is 
easy to verify. Therefore two of its three eigenvalues are equal to zero, while the remaining 
non-zero eigenvalue represents the solution we are looking for. 

The analytic expression of the third eigenvalue can be obtained by using the property 
that the trace of the 3x3 matrix C^^ ■ A is equal to the sum of its three eigenvalues, and 
in our case to the eigenvalue we are looking for. From these considerations we derive the 
following expression for the optimized signal-to-noise ratio SNR^^^^ : 

SNRLp,= r 4'^* {C-% df . (71) 



We can summarize the results derived in this section, which are given by equations (jHEEJ), 
in the following way: 

(i) among all possible interferometric combinations LISA will be able to synthesize with its 
four generators a,/5,7,C, the particular combination giving maximum signal-to-noise ratio 
can be obtained by using only three of them, namely 7); 

(ii) the expression of the optimal signal-to-noise ratio given by equation (ffT|l implies that 
LISA should be regarded as a network of three interferometer detectors of gravitational 
radiation (of responses (a,/3, 7)) working in coincidence jiol. IsS^. 



A. General application 

As an application of equation (ffT|) . here we calculate the sensitivity that LISA can reach 
when observing sinusoidal signals uniformly distributed on the celestial sphere and of random 
polarization. In order to calculate the optimal signal-to-noise ratio we will also need to use a 
specific expression for the noise correlation matrix C. As a simplification, we will assume the 
LISA arm-lengths to be equal to its nominal value L = 16.67 sec, the optical-path noises to 
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be equal and uncorrelated to each other, and finally the noises due to the proof-mass noises 
to be also equal, uncorrelated to each other and to the optical-path noises. Under these 
assumptions the correlation matrix becomes real, its three diagonal elements are equal, and 
all the off-diagonal terms are equal to each other, as it is easy to verify by direct calculation 
[isj ]. The noise correlation matrix C is therefore uniquely identified by two real functions, 
Sa and Sai3, in the following way 



C 



^ Sa Sap Sap ^ 

Sa(3 Sa Sap 
y Sap Sap Sa j 



The expression of the optimal signal-to-noise ratio assumes a rather simple form if we 
diagonalize this correlation matrix by properly "choosing a new basis" . There exists an or- 
thogonal transformation of the generators (5, /5, 7) which will transform the optimal signal- 
to-noise ratio into the sum of the signal-to-noise ratios of the "transformed" three inter- 
ferometric combinations. The expressions of the three eigenvalues (which are real) 

of the noise correlation matrix C can easily be found by using the algebraic manipulator 
Mathematica [2^, and they are equal to 

(J'l = 1^2 = Sa — Sap , fl^ = Sa + 2 Sap ■ (72) 

Note that two of the three real eigenvalues, (/ii,/i2), are equal. This implies that the eigen- 
vector associated to /is is orthogonal to the two-dimensional space generated by the eigen- 
value /ii, while any chosen pair of eigenvectors corresponding to /ii will not necessarily be 
orthogonal. This inconvenience can be avoided by choosing an arbitrary set of vectors in 
this two-dimensional space, and by ortho-normalizing them. After some simple algebra, we 
have derived the following three ortho-normalized eigenvectors 

vi = ^ (-1,0,1) , V2 = -^ (1,-2,1) , V3 = -^ (1,1,1). (73) 

Equation ()73|) implies the following three linear combinations of the generators (5, (3, 7) 

A=-^(7-5) , E = -^ (5-2^+7) , T=i= (5 + ^ + 7), (74) 

where A, E, and T are italicized to indicate that these are "orthogonal modes" . Although the 
expressions for the modes A and E depend on our particular choice for the two eigenvectors 
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(vi,V2), it is clear from our earlier considerations that the value of the optimal signal-to- 
noise ratio is unaffected by such a choice. From equation ()74|) it is also easy to verify that 
the noise correlation matrix of these three combinations is diagonal, and that its non-zero 
elements are indeed equal to the eigenvalues given in equation (j72j) . 

In order to calculate the sensitivity corresponding to the expression of the optimal signal- 
to- noise ratio, we have proceeded similarly to what was done in 11 1, 2, and described 
in more detail in 3^. We assume an equal-arm LISA {L = 16.67 light seconds), and 
take the one-sided spectra of proof mass and aggregate optical-path-noises (on a single 
link), expressed as fractional frequency fluctuation spectra, to be 0); 8^'^°°^ = 

2.5 X 10-^^ [f/lHzY^ Hz"^ and p'^*^ = 1.8 x lO^^^ [f I^Hzf Hz"\ respectively. We 

also assume that aggregate optical path noise has the same transfer function as shot noise. 

The optimum SNR is the square root of the sum of the squares of the SNRs of the three 
"orthogonal modes" {A, E, T). To compare with previous sensitivity curves of a single LISA 
Michelson interferometer, we construct the SNRs as a function of Fourier frequency for 
sinusoidal waves from sources uniformly distributed on the celestial sphere. To produce the 
SNR of each of the {A, E, T) modes we need the gravitational wave response and the noise 
response as a function of Fourier frequency. We build up the gravitational wave responses 
of the three modes {A,E,T) from the gravitational wave responses of (a,/3,7). For 7000 
Fourier frequencies in the ~10~^ Hz to ~1 Hz LISA band, we produce the Fourier transforms 
of the gravitational wave response of (a,/3,7) from the formulas in 11 1, j^^. The averaging 
over source directions (uniformly distributed on the celestial sphere) and polarization states 
(uniformly distributed on the Poincare sphere) is performed via a Monte Carlo method. 
From the Fourier transforms of the (a,/?, 7) responses at each frequency, we construct the 
Fourier transforms of (A, E, T). We then square and average to compute the mean-squared 
responses of {A, E, T) at that frequency from 10^ realizations of (source position, polarization 
state) pairs. 

We adopt the following terminology: we refer to a single element of the module as a data 
combination; while a function of the elements of the module, such as taking the maximum 
over several data combinations in the module or squaring and adding data combinations 
belonging to the module, is called as an observable. The important point to note is that 
the laser frequency noise is also suppressed for the observable although it may not be an 
element of the module. 
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The noise spectra of {A, E, T) are determined from the raw spectra of proof-mass and 
optical-path noises, and the transfer functions of these noises to {A, E, T). Using the transfer 
functions given in jl3|, the resulting spectra are equal to, 

SAif) = SEif) = 16 sin2(vr/L) [3 + 2cos(27r/L) + cos(47r/L)]^^P™''^ --^(/) 

+ 8 sin^n fL) [2 + cos(27r/L)] S^*'™' ^^'^f) , (75) 



Srif) = 2[l + 2cos(27r/L)]2 [4 sm\7rfL)SP"'"f + S^p''""^ f)] . (76) 

Let the amplitude of the sinusoidal gravitational wave be h. The SNR for, e.g. A, SNR^i, at 
each frequency / is equal to h times the ratio of the root-mean-squared gravitational wave 
response at that frequency divided by a/5'a(/) B, where B is the bandwidth conventionally 
taken to be equal to 1 cycle per year. Finally, if we take the reciprocal of SNR^i//;, and 
multiply it by 5 to get the conventional SNR = 5 sensitivity criterion, we obtain the sen- 
sitivity curve for this combination which can then be compared against the corresponding 
sensitivity curve for the equal-arm Michelson Interferometer. 

In Figure El we show the sensitivity curve for the LISA equal-arm Michelson response 
(SNR = 5) as a function of the Fourier frequency, and the sensitivity curve from the opti- 
mum weighting of the data described above: 5h/ a/SNR^ + SNR|; + SNR^. The SNRs were 
computed for a bandwidth of 1 cycle/year. Note that at frequencies where the LISA Michel- 
son combination has best sensitivity, the improvement in signal-to-noise ratio provided by 
the optimal observable is slightly larger than 

In Figure m we plot the ratio between the optimal SNR and the SNR of a single Michelson 
interferometer. In the long- wavelength limit, the SNR improvement is V2. For Fourier fre- 
quencies greater than or about equal to 1/L, the SNR improvement is larger and varies with 
the frequency, showing an average value of about y/S. In particular, for bands of frequencies 
centered on integer multiples of 1/L, SNRt contributes strongly and the aggregate SNR in 
these bands can be greater than 2. 

In order to better understand the contribution from the three different combinations to 
the optimal combination of the three generators, in Figure |H1 we plot the signal-to-noise 
ratios of {A,E,T) as well as the optimal signal-to-noise ratio. For an assumed h = 10^^^, 
the SNRs of the three modes are plotted versus frequency. For the equal-arm case computed 
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FIG. 6: The LISA Michelson sensitivity curve (SNR = 5) and the sensitivity curve for the optimal 
combination of the data, both as a function of Fourier frequency. The integration time is equal to 
one year, and LISA is assumed to have a nominal armlength L = 16.67sec. 




log.,Q(frequency, Hz) 

FIG. 7: The optimal SNR divided by the SNR of a single Michelson interferometer, as a function 
of the Fourier frequency /. The sensitivity gain in the low- frequency band is equal to \/2, while it 
can get larger than 2 at selected frequencies in the high-frequency region of the accessible band. 
The integration time has been assumed to be one year, and the proof mass and optical path noise 
spectra are the nominal ones. See the main body of the paper for a quantitative discussion of this 
point. 



here, the SNRs of A and E are equal across the band. In the long wavelength region of the 
band, modes A and E have SNRs much greater than mode T, where its contribution to the 
total SNR is negligible. At higher frequencies, however, the T combination has SNR greater 
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FIG. 8: The SNRs of the three combinations, {A, E, T), and their sum as a function of the Fourier 
frequency /. The SNRs of A and E are equal over the entire frequency band. The SNR of T is 
significantly smaller than the other two in the low part of the frequency band, while is comparable 
to (and at times larger than) the SNR of the other two in the high-frequency region. See text for 
a complete discussion. 

than or comparable to the other modes and can dominate the SNR improvement at selected 
frequencies. Some of these results have also been obtained in [igj ]. 



B. Optimization of SNR for binaries with known direction but with unknown 
orientation of the orbital plane 

Binaries will be important sources for LISA and therefore the analysis of such sources is of 
major importance. One such class is of massive or super massive binaries whose individual 
masses could range from IO^Mq to 10^ Mq and which could be upto a few Gpc away. 
Another class of interest are known binaries within our own galaxy whose individual masses 
are of the order of a solar mass but are just at a distance of a few kpc or less. Here the 
focus will be on this latter class of binaries. It is assumed that the direction of the source 
is known, which is so for known binaries in our galaxy. However, even for such binaries, 
the inclination angle of the plane of the orbit of the binary is either poorly estimated or 
unknown. The optimization problem is now posed differently: the SNR is optimized after 
averaging over the polarizations of the binary signals, so the results obtained are optimal 
on the average, that is, the source is tracked with an observable which is optimal on the 
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average For computing the average, a uniform distribution for the direction of the 

orbital angular momentum of the binary is assumed. 

When the binary masses are of the order of a solar mass and the signal typically has a 
frequency of a few mHz, the GW frequency of the binary may be taken to be constant over 
the period of observation, which is typically taken to be of the order of an year. A complete 
calculation of the signal matrix and the optimization procedure of SNR, is given in [2^. 
Here we briefly mention the main points and the final results. 

A source fixed in the Solar System Barycentric reference frame in the direction {0b,4>b) 
is considered. But as the LISA constellation moves along its heliocentric orbit, the apparent 
direction {9i,(f)L) of the source in the LISA reference frame (xr,, V t,. zj) changes with time. 
The LISA reference frame {x^^ yi, zi) has been defined in reference j20| as follows: the origin 
lies at the center of the LISA triangle and the plane of LISA coincides with the [xl-, Ul) plane 
with spacecraft 2 lying on the xl axis. Fig. (jH)) displays this apparent motion for a source 
lying in the ecliptic plane, that is with 9b = 90° and 0^ = 0°. The source in the LISA 
reference frame describes a figure of 8. Optimizing the SNR amounts to tracking the source 
with an optimal observable as the source apparently moves in the LISA reference frame . 

Since an average has been taken over the orientation of the orbital plane of the binary or 
equivalently over the polarizations, the signal matrix A is now of rank 2 instead of rank 1 
as compared with the application in the previous subsection. The mutually orthogonal data 
combinations A, E, T are convenient in carrying out the computations because in this case 
as well, they simultaneously diagonalize the signal and the noise covariance matrix. The 
optimization problem now reduces to an eigenvalue problem with the eigenvalues being the 
squares of the SNRs. There are two eigen- vectors which are labelled as -u+^x belonging to 
two non-zero eigenvalues. The two SNRs are labelled as, SNR+ and SNRx corresponding 
to the two orthogonal (thus statistically independent) eigen- vectors w+^x- As was done in 
the previous subsection the two SNRs can be squared and added to yield a network SNR, 
which is defined through the equation: 

SNRLwork = SNR^ + SNRl . (77) 

The corresponding observable is called as the network observable. The third eigenvalue is 
zero and the corresponding eigenvector orthogonal to and gives zero signal. 

The eigenvectors and the SNRs are functions of the apparent source direction parameters 
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FIG. 9: Apparent position of the source in the sky as seen from LISA frame for 
{Ob = 90°, (f)B = 0°). The track of the source for a period of one year is shown on the unit 
sphere in the LISA reference frame. 
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4'l) in the LISA reference frame, which in turn are functions of time. The eigenvectors 
optimally track the source as it moves in the LISA reference frame. Assuming an observation 
period of an year, the SNRs are integrated over this period of time. The sensitivities are 
computed according to the procedure described in the previous subsection IVI Al The results 
of these findings are displayed in FigfTUl It shows the sensitivity curves of the following 
observables: 

1. The Michelson combination X (faint solid curve). 

2. The observable obtained by taking the maximum sensitivity among X, Y and Z for 
each direction, where Y and Z are the Michelson observables corresponding to the 



remaining two pairs of arms of LISA This maximum is denoted by max[X, Y, Z] 
(dash-dotted curve) and is operationally given by switching the combinations X, Y, Z 
so that the best sensitivity is achieved. 

3. The eigen-combination {7+ which has the best sensitivity among all data combinations 
(dashed curve). 
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FIG. 10: Sensitivity curves for the observables: Michelson, max[X, y, Z], and network for the 
source direction {6b = 90°, 0b = 0°). 
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4. The network observable (solid curve). 

It is observed that the sensitivity over the band-width of LISA increases as one goes 
from (1) to (4). Also it is seen that the max[X, F, Z] does not do much better than X. 
This is because for the source direction chosen 6b = 90°, X is reasonably well oriented and 
switching to Y and Z combinations does not improve the sensitivity significantly. However, 
the network and v+ observables show significant improvement in sensitivity over both X 
and max[X, Y, Z]. This is the typical behavior and the sensitivity curves (except X) do not 
show much variations for other source directions and the plots are similar. Also it may be 
fair to compare the optimal sensitivities with max[X, Y, Z] rather than X. This comparison 
of sensitivities is shown in Fig lTTl where the network and the eigen-combinations -u+^x are 
compared with max[X, Y, Z\. 

Defining: 

SNR,(/) 

<^a{f) = ^ 7TT , (78) 

where the subscript a stands for network or +, x and SNRinax[x,y,z] is the SNR of the 
observable max[X, y, Z], the ratios of sensitivities are plotted over the LISA band-width. 
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FIG. 11: Ratios of the sensitivities of the observables network, v+^x with max[X, y, Z] for the 
source direction 63 = 90°, 0b = 0°. 
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The improvement in sensitivity for the network observable is about 34% at low frequencies 
and rises to nearly 90% at about 20 mHz, while at the same time the -0+ combination shows 
improvement of 12% at low frequencies rising to over 50% at about 20 mHz. 



VII. CONCLUDING REMARKS 



In this article we have summarized the use of TDI for canceling the laser phase noise from 
heterodyne phase measurements performed by a constellation of three spacecraft tracking 
each other along arms of unequal length. Underlying the TDI technique is the mathematical 
structure of the theory of Grobner basis and the algebra of modules over polynomial rings. 
These methods have been motivated and illustrated with the simple example of an unequal 
arm interferometer in order to give a physical insight of TDI. Here, these methods have been 
rigorously applied to the idealized case of a stationary LISA for deriving the generators of 
the module from which the entire TDI data set can be obtained; they can be extended in a 
straightforward way to more than three spacecraft for possible LISA follow-on missions. The 
stationary LISA case was used as a propaedeutical introduction to the physical motivation of 
TDI, and for further extending it to the realistic LISA configuration of free-faUing spacecraft 
orbiting around the Sun. The TDI data combinations cancehng laser phase noise in this 
general case are referred to as second generation TDI, and they contain twice as many terms 
as their corresponding first generation combinations valid for the stationary configuration. 

As a data analysis application we have shown that it is possible to identify specific TDI 
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combinations that will allow LISA to achieve optimal sensitivity to gravitational radiation 



m 



19|,|20|. The resulting improvement in sensitivity over that of an unequal-arm Michelson 
Interferometer, in the case of monochromatic signals randomly distributed over the celestial 
sphere and of random polarization, is non negligible. We have found this to be equal to 
a factor of in the low-part of the frequency band, and slightly more than in the 
high-part of the LISA band. The SNR for binaries whose location in the sky is known but 
their polarization is not can also be optimized, and the degree of improvement depends on 
the location of the source in the sky. 

As a final remark we would like to emphasize that this field of research, TDI, is still very 
young and evolving. Possible physical phenomena, yet unrecognized, might turn out to be 
important to account for within the TDI framework. The purpose of this review was to 
provide the basic mathematical tools needed for working on future TDI projects. We hope 
to have accomplished this goal, and that others will be stimulated to work in this new and 
fascinating field of research. 



Acknowledgement 

SVD acknowledges support from IFCPAR, Delhi, India under which the work was carried 
out in collaboration with J-Y Vinet. This research was performed at the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with the National Aeronautics 
and Space Administration. 

APPENDIX A: GENERATORS OF THE MODULE OF SYZYGIES 

We require the 4-tuple solutions (gs, q[, gj, 13) to the equation: 

(1 - xyz)qs + {xz - y)q[ + x{l - z^y^ + (1 - x')?^ = , (Al) 

where for convenience we have substituted x = Vi,y = V2, z = V^. q3,q'i,q2,q'3 are polyno- 
mials in x,y,z with integral coefficients i.e. in Z[x,y,z]. 
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We now follow the procedure in the book by Becker et al. 
Consider the ideal in Z[x,y,z] (or Q[x,y,z] where Q denotes the field of rational num- 
bers) , formed by taking linear combinations of the coefficients in Eq. ()Alj) /i = 1 — xyz, f2 = 
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xz — y, h — x{l — z^), /4 = 1 — x'^. A Grobner basis for this ideal is: 

^ = {s-i = - 1, 5-2 = - 1, 5-3 = a; - yz] . 



(A2) 



The above Grobner basis is obtained using the function GroebnerBasis in Mathematica. One 
can check that both the = 1, 2,3,4 and gj,j = 1,2,3 generate the same ideal because 
we can express one generating set in terms of the other and vice- versa: 



9j ~ Cjifi , 



(A3) 



where d and c are 4x3 and 3x4 polynomial matrices respectively, and are given by, 



/-I 



d 



\ 



^ -X z'^-1^ 



\ 



-1 -y 
z 1 



(A4) 



-yz 
y z 
-X 
\-l -(x + yz)J 

The generators of the 4-tuple module are given by the set A \^B* where A and B* are the 
sets described below: 

A is the set of row vectors of the matrix I — d-c where the dot denotes the matrix product 
and I is the identity matrix, 4 x 4 in our case. Thus, 

Qi — {z^ — 1, 0, X — yz, 1 — z"^) , 

as = {0,z{l- z'^),xy - z,y{l - z^)) , 

aa = {0,0,1 - x\x{z^ -1)), 

a4 = {—z'^,xz,yz,z'^) 



(A5) 



We thus first get 4 generators. The additional generators are obtained by computing the S- 
polynomials of the Grobner basis Q. The S-polynomial of two polynomials gi, g2 is obtained 
by multiplying gi and g2 by suitable terms and then adding, so that the highest terms cancel. 
For example in our case gi — z"^ — 1 and g2 — y^ — 1 and the highest terms are z'^ for gi and 
y'^ for g2 . Multiply gi by y"^ and g2 by z'^ and subtract. Thus, the S-polynomial pu of gi 
and g2 is: 

2 2 2 2 ^ \ c\ 

Pi2 = y 91- z g2 = z -y . (A6) 

Note that order is defined {x » y » z) and the y'^z'^ term cancels. For the Grobner 
basis of 3 elements we get 3 S-polynomials Pi2,Pi3,P23- The Pij must now be re-expressed in 
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terms of the Grobner basis Q. This gives a 3 x 3 matrix b. The final step is to transform to 
4-tuples by multiplying b by the matrix c to obtain b* — b ■ c. The row vectors b^^i— 1, 2, 3 
of b* form the set B*: 

bl = {z'-l,y{z'-l),x{l-y'),(y'-i)(z'-l)), 

b; = {0,z{l-z^),l-z'-x{x-yz),{x-yz){z'-l)), 

K = {-x + yz,z-xy,l-y^,0). ( A7) 

Thus we obtain 3 more generators which gives us a total of 7 generators of the required 
module of syzygies. 

APPENDIX B: CONVERSION BETWEEN GENERATING SETS 

We list the three sets of generators and relations among them. We first list below a, /3, 7, (: 

a = (-1, -z, -xz, l,xy,y) , 

P = {-xy, -1, -x,z, l,yz) , 

7 = {-y,-yz,-i,xz,x,i) , 

C = (-x,-y,-z,x,y,z) . (Bl) 
We now express the and bj in terms of a, 7, C: 

«i = 7-<, 
a2 — a — z(3 , 

as — —za + /3 — x^i + xzC, , 

04 = zC, 

b\ = -ya + + 7 - <, 

b\ = {l-z^)p-xj + xzC, 

bl^P-yC. (B2) 
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Further we also list below a, /3, 7, C in terms of X^-^'>: 

a = , 
(5 = X(^) , 

C = . (B3) 

This proves that since the ai^h* generate the required module, the q;,/?,7, C and 
X'^^\ A = 1, 2, 3, 4 also generate the same module. 

The Grobner basis is given in terms of the above generators as follows: 
= C, = X('\ (5(3) = ^(4) = a and G^^) = 03 
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